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Abstract 

The behavior of the most probable values of the order parameter x and the amplitude (j> of conjugate 
force fluctuations is studied for a stochastic system with a colored multiplicative noise with absorbing states. 
The phase diagrams introduced as dependencies the noise self-correlation time vs temperature and noise 
growth velocity are defined. It is shown that phase half-plane (a;, <j)) can be split into isolated domains of 
large, intermediate, and small values of x. System behavior in these domains is studied by the probability 
represented as path integral. In the region a; ^ 1, the trajectories converge to the point x = = for 
< a < 1/2 and to s = 0, oo for 1/2 < a < 1. In the former case, the probability of realization of 

trajectories is finite, while in the latter case it is vanishingly small, and an absorbing state can be formed. 

A wide class of stochastic systems with multiplicative noise whose intensity is a function of the stochastic 
variable x determining the state of the subsystem can exhibit absorbing states. An absorbing configuration is 
one in which a system can get trapped, from which it can not escape The most popular example of this kind 
is the Verchiilst model proposed for explaining the population kinetics Q]; the analysis of this model remains a 
topical question Q . Another example is the percolation model Q in which a linear dependence on x is observed 
for noise intensity. Indicated models can be generalized by representing the noise intensity as a stochastic 
variable power dependence with an arbitrary exponent. A kinetic features of stochastic system with that kind 
noise was studied in |5| in the white noise approximation and with accounting inhomogeneity of the space. 
An evolution of averages in that model was discussed in Q. A distinguishing feature of our approach is that 
we shall study the dependencies of most probable values vs system parameters using supposition of the noise 
correlation. This is realized on the basis of a field scheme in which the probability functional has an exponential 
form with the exponent that can be reduced (correct to the sign) to the standard action 0. Consequently, a 
description of phase transition in terms of the most probable values corresponds to the application of the least 
action principle. Such an approach makes it possible to extend the self-consistent field method to a description 
of systems with multiplicative noise. This reveals the following nontrivial peculiarity: the most probable values 
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can tend either to free energy minima or to points corresponding to unstable thermodynamic states depending 
on the system parameters. The singular nature of noise indicating the tendency of its intensity to zero for the 
value x = of the order parameter (more probably magnitude of x) is responsible for the emergence of a region 
on the phase portrait in which all trajectories converge to the axis x — over a finite time. In other words, an 
absorbing state is formed, in which the system becomes closed 

We consider following model. The time dependence of the hydrodynamic mode amplitude x(t) is determined 
by the Langevin equation 

x = fo{x)+go{x)C{t). (1) 

Here the dot indicates the differentiation with respect to time t, deterministic force /o(a;) is defined according to 
the thermodynamic potential Vo{x) as fa = -~dVo{x) / dx^ C{t) is the colored multiplicative noise with amplitude 
defined by function 50(2;). A specific form for the noise we choose in kind of Ornstein-Ulenbeck process, 
i.e. 

rC = -C + ^W, (2) 

where r is a self-correlation time, C is the Gaussian distributed variable with zero mean and exponentially 
decaying correlations: 

(CWC(i'))-(l/2r)exp(-|t-t'|/r), (3) 

^{t) is a white noise — namely, Gaussian stochastic variable with zero mean and ^-correlated: {^{t)^{t')) = 
S(t — t'). The amplitude go{x) of the multiplicative noise in Eq.(|I]) is chosen in the simple form 

9i^) = kr, (4) 

where a is a positive exponent varying from to 1 We shall describe our system by master parameter, 
self-correlation time and velocity of noise growth defined by exponent a. 

If we take time derivative of Eq.(^, replace first C in terms of C and ^ from Eq.(H) and then C in terms of x 
and x from Eq.(|l|), we obtain the following non-Markovian stochastic differential equation (SDE) (sec Rcf. |^): 

r (x - x'^dxgo{x)/go{x)) = -a{x)x + /o(.t) + go{x)^{t), (5) 

where 



a{x) = 1 - Tfo{x)dxh-i 



fo{x) 



(6) 



(throughout the paper, the Ito interpretation for the SDE will be meant). In the limit = a system behavior 
was described in |l^, |ll|. Here, following "unified colored noise approximation" Ref. we can recover a 
Markovian SDE. It needs to use adiabatic elimination (neglecting x) and to neglect so that the system's 
dynamics be governed by a Fokker-Planck equation. In order to make correct transition to an ordinary, linear 
in X, SDE we have to use Ito's differential rule |l3|. In above mentioned suppositions SDE (5) takes form 

a{x)x = fo{x)+g„{x)at). (7) 
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Expression in Ihs in Eq.(|^) defines time derivative for variable z. The relation between z, and x is set by 
expression dz — a{x)dx. An equivalent equation to Eq.(^ can be write down as 

f(x) + h{x)+g{x)m. (8) 

with 



m = Hx) --,[^) a.ln.(.), = ^. (9) 

Here we stress that drift term h{x) appears only by accounting the noise self-correlation time. 
Let us now go over to a new field y(t) connected with the initial field x{t) through the relation 

^=9{x{y)). (10) 

The new stochastic field taking into account Eq.(|l0|) satisfies the Langevin equation 

y = h{x{y))+^, (11) 

with additive noise and effective force 

h = {f + h)/g-d^g/2. (12) 

The obtained equation (|l^) allows us to use the standard field scheme [Q based on analysis of the generating 
functional. The latter has the form of the functional Laplace transform 

Z{u{t)} = J Z{y{t)}exp (^J uyd?j Vy{t) (13) 

for the partition function 



(14) 



Z{y{t)} = {^5[y-h~i]dei - 

where Vy denotes integration over all paths starting at j/(0) for t = Q and ending at y{tf) for t = tf. It is 
defined as 

N-l 

Vy=\im l[dy{U), (15) 

wlu'vv y{ti) is a field at time ti = ie, having sliced the interval to tf in N parts of size e — tf/N. The 
argument of the (5- function in Eq.(|l^) can be reduced to the Langevin equation (pj]), and the determinant 
ensuring a transition from continual integration with respect to ^{t) to y{t) is equal to unity in Ito calculus. 

In the approach developed by Zinn-Justin n-fold variation of functional (^|) over the auxiliary field u{t) 
makes it possible to find the n-th order correlator for the hydrodynamic mode amplitude y{t) and to construct 
the perturbation theory. We shall proceed, however, from expression ( |l^ ) for the conjugate functional Z{y{t)} 
whose variation leads to the most probable realization of the stochastic field y{t). 

Going over to an analysis of functional (p^), we write the 5- function in integral form 



S{y{t)}^ / exp - / qydt]Vq. (16) 
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Averaging over the noise with the help of the Gaussian distribution 



(17) 



corresponding to (g) and taking into account Eq.(|l6|), we reduce functional Eq.(14) to the standard form 

Z{y{t)}^ J P{y{t),q{t)}Vq, P = e-^. (18) 
Here the probability distribution P{y, q} is specified by the action S = J Cdt, where the Lagrangian is given by 

C{y,q) = q(y-h) -qy2. (19) 



We shall use in subsequent analysis the Euler equations 

dC _ d_dC _ m 
dz dt dz dz ' 

in which the dissipative function is defined as 



(20) 



(21) 



As a result, the equations for the most probable realizations of the stochastic fields y{t) and q{t) assume the 
form 



y^h + q, 
q = —q\l + dxh ) — h. 



(22) 
(23) 



A comparison of Eq.(^2|) with the stochastic equation (11) having the same form shows that the fields y{t) 
and q{t) are the most probable values of amplitudes of the auxiliary hydrodynamic mode defined by relation 
(^ as well as by fluctuation of the conjugate force. Obviously, the latter can be reduced to the conjugate 
momentum q = dC/dy. 

In order to return to the initial stochastic field x{t), we shall use the relation ( p^ and definition (p^). As a 
result, the Lagrangian (19) and the dissipative function ( pl| ) assume the form 



n{x) = xy2g\ 

where we have used the definition of the conjugate momentum = dC/dx leading to the relation 

= q/g- 

In this case, the Euler equations ( ^ ) assume the form 

x^f + h- gdxg/2 + g^ct), 
= -^[l + d^f + d,h - dx{gdxg)/2 + ^gd^g] - (/ + h)/g'' + dxg/2g. 



(24) 
(25) 

(26) 

(27) 
(28) 
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Similar equations can be obtained directly from the system (p2|), ( p3| ) using the relations and (|2^) as well 
as definition 

In order to calculate noise self-correlation time influence on the system behavior we shall use the x^-model 
for potential 

Vr,{x) = '-x' + \x\ e e [-1, 1] (29) 

and the definition (^ of multiplicative function. Since the Lagrangian ( p4|) does not change its form upon 
simultaneous reversal of the signs of x and the phase portraits will possess central symmetry relative to the 
origin a; = = 0. On the other hand, the axis a; = on which the noise intensity assumes zero value is singular 
in accordance with (^) . For this reason, we can confine our analysis only to the upper part of the phase plane 
corresponding to the value of a; > 0. 

First of all let us investigate steady states of the system. In the stationary case i = = 0, we have two 
equations 



' 9l +2 



2 i r,3xlnffo, (30) 



l^^j^ln 


'fo 


I % 


.9o_ 



+ hd,\na)d,lng„~^}^0. (31) 



From Eq.(|30|) it follows that at noise self-correlation time not exceeding the value tq (Fig.l) the form of the 
phase portrait is characterized by the presence of a single saddle point S whose position is defined by the 
solution of Eq.(^Il) and corresponds to cf) ^ 0. Figure 2 shows that bifurcation takes place at = at noise 
self-correlation time tq (Fig. 2a, 2b). This is accompanied by the emergence of an additional saddle point S and 
an attractive node C whose positions are determined by the condition ip — and specified by the coordinates 
x^. 

In Fig. 3 we show critical magnitude of the noise self-correlation time Tc vs temperature e (as maximum of 
the dependence To(a)). As shown in Fig.l, above Tc the system is always ordered for any velocity of the noise 
growth, defined by exponent a. Obviously, the saddle point S and the node C merge at the point corresponding 
to the correlation time Tc- 

The above analysis shows that in systems with r < tq the stationary state corresponds to the point S whose 
coordinates are defined as 

<^ = -^ + Ja.5, (32) 

g 2 

dxif + h) = -(l)gd.j;g + ^d^igd^g]. (33) 

The meaning of these coordinates becomes clear if we proceed to the additive limit g{x) I: condition ( ^^ 
indicates that the most probable value ((> of the fluctuation amplitude of conjugate force for homogeneous systems 
has the sign opposite to that of / + according to Eq.(^), the "susceptibility" (d'^V)^^ = —l/dx{f + h) in this 
case assumes an infinitely large value. Thus, node S corresponds to the stationary state of a thermodynamic 
system which is unstable to a transition to the ordered phase. Then the phase portrait has the form shown in 
Fig.4a. 
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At temperatures t > tq, the steady state with the coordinates 



^ = 0, 



(34) 



2cr/o = godx9o- 



(35) 



is formed at the point C. In the additive hmit, this point corresponds to the state of thermodynamic equihbrium. 
Figure 4 shows that the corresponding phase portrait can be obtained not only upon an decrease in temperature 
(see Fig. 4a). Here we stress that the system can pass from equihbrium (ordered) domain trough unstable 
(disordered) to stable if we increase the velocity of the noise growth (see Fig.l and Fig. 4b). 

Naturally, a real thermodynamic system in the process of its evolution tends to the equilibrium state Eqs.(|3^, 
^5|) rather than to the unstable state Eqs.(|3|, |3^). The equilibrium state corresponds to small and large values 
of the exponent a for any values of the parameter t. The domain of medium values of the exponent a at small 
noise correlations corresponds to unstable states. 

Comparing Fig. 4a, Fig. 4b we see that they are distinguished by the location of the attraction node on the axis 
= 0. The general form of phase portraits shown in Fig. 4 is characterized by the presence of two separatrices 
with branches PQ and MSqN. They divide the phase plane into three isolated regions corresponding to large, 
intermediate, and small values of the order parameter x. The first region is characterized by an indefinite 
increase in the values of the quantities x and with time t ^ oo. It will be sown below that this is not 
realized in actual practice. The region of intermediate values of x in which the system ultimately goes over to a 
stationary ordered state is most interesting. It is this region that determines the phase transition kinetics. The 
formation of the region corresponding to values x ^ 1 is associated with the multiplicative nature of noise. In 
this region, the order parameter x{t) tends with time to the value x = 0. 

Let us analyze the behavior of the system in each of these regions. For this purpose, we consider the 
probability of realization of a phase trajectory corresponding to different initial values xq = x{t = 0). In 
accordance with (|l^), the probability can be written in the form 



where expressions ( p4|) and (g7|) are taken into account and integration is carried out along the corresponding 
trajectory. The dependence P{xq) obtained for an exponent a < 1/2 is shown in Fig. 5 (curves 1, 2). Apart from 
the trivial increase in probability ( |3^ ) which approaches the origin for values of xq corresponding to separatrices, 
the jumps near which the value of P{xo) can increase insignificantly are observed. Outside the region bounded 
by the (outer) separatrix, we have P = since x(t), <j)[t) ±oo for i cx) in this case. 

Such a behavior of the probability P{xq) can be explained by the form of the time dependencies x{t) and 
4>{t) during relaxation of the initial value. Far away from the region of 5*0 (see Fig. 4a), the quantities and 
X rapidly change their values, the change slowing down as we approach this region. Such a behavior can be 
explained by the fact that the action S{x{t), (/'(t)} changes much more slowly near this region than away from 
it. This can be visualized by associating the region SqCS with the bed of a large river |15[|. 




(36) 
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For values of the exponent a > 1/2, the integrand in formula (|36|) diverges, and the probability P assumes 
zero value for «C 1. Typically, this divergence is observed only in the region of phase portrait bounded by 
the separatrix branch SoO in Fig. 4b. 

In order to explain the form of the dependence P{xo), we analyze the behavior of the quantities x{t) and 
(f>{t) for various values of exponent a. For this purpose, we put (p = in Eq.(|2^). The obtained quadratic 
equation gives stationary values of conjugate momentum in the limit a; ^ 0: 

(37) 

Thus, for a < 1/2, the system tends with time to the origin x — (f> — 0, and the attraction node jumps to 
infinity {(f) — > oo, x = 0) as a exceeds the critical value a = 1/2. The corresponding integrand in distribution 
il),i.e., 




is characterized by the sign inversion in the exponent upon a transition through the critical value a = 1/2. 
Substituting Eq.(^^ for a < 1/2 into the equation of motion ( p7| ) and retaining the leading term in it, we 
obtain the equation 2x — —ax'^"'^^ which gives the time dependence of the order parameter: 

= a(l -a)(fo -<), t<tQ, a<l/2, (39) 

where to is the integration constant defining the time during which the point gets to the axis x = 0. The 
substitution of dependence Eq. (|39|) into Eq.(|3^) at a < 1/2 and of the obtained expression into the integral in 
Eq.(^6|) shows that the probability P{xq) of realization of a trajectory in the region x -C 1 of the phase portrait 
differs from zero (see Fig. 5, curves 1, 2). 

A completely different situation is observed for an exponent a > 1/2. In this case, Eq.([27|) can be reduced 



to 2x = —(1 — a)x'^"'^^, leading to Eq.(39) in which the factor a(l — a) is replaced by (1 
for a > 1/2 the expression g'^cjP acquires an exponent with the opposite sign in accordance with Eq.(p8|) so 
that it assumes an infinitely large value for < — > io- As a result, the probability Eq.(|3^) becomes vanishingly 
small (see Fig. 5, curves 3, 4). The physical reason behind such a behavior is that the system gets to the axis 
a; = over a finite time interval to < oo, which ensures an infinitely large value of the conjugate momentum 
(j) (X x~^ oc (to — t)~^^^^^~°'\ This can be visualized as the precipitation of condensate (absorbing state) of 
configuration points from the phase portrait domain x <^ 1 onto the abscissa axis for — > oo. Note that the 
condition to < cxo is satisfied only below the separatrix branch ^oO in Fig. 4b, while in the region NSqO we 
have to = oo, and the convergence of the integrand in (^) is not manifested. Consequently, the equality P = 
holds only below the curve SqO (see Fig. 5 and Fig. 4b). 

In accordance with (|39|), the time dependence of the most probable magnitude x{t) oc is defined by the 
exponent H^^ = 2(1 — a) whose magnitude determines the fractional dimension D = characterizing the 
domain a; ^ 1 of the phase portrait of a system with multiplicative noise pM. In the additive limit a — 0, we 
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have dimension D = 2oi the phase plane as expected. This means that as the time t ^ oo, the phase trajectories 
of the system fill the entire phase plane. The increase of the exponent a > leads to a decrease in the fractional 
dimension D which assumes the critical value £) = 1 for a = 1/2. As the value of a increases further, the 
fractional dimension of the set of points on the plane x, t, which is the law of motion x{t), becomes smaller than 
unity. The physical reason behind such a behavior is the above mentioned absorption of configuration points 
by the axis a; = for — > oo. 

Finally, it must be stressed that the noise correlation time can play role of master parameter as the temper- 
ature and produces transition to the ordered state. Moreover we shown that the increasing of the velocity of the 
noise intensity growth produces appearance of the domain of disordered states. This domain disappears if noise 
self-correlation time does not exceed value Tc- Path integral solution shown that absorbing state appearance is 
characterized by the vanishingly small probability of the system states realization. Asymptotic time dependence 
of the order parameter x explains picture of absorbing states by the fractal dimension smaller than unity. 
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FIGURE CAPTIONS 



Fig.l. Phase diagram indicating appearance of the ordered phase (x ^ 0); noise correlation time r vs 
exponent a (curves 1, 2 correspond to e = 0.65, 0.7). 

Fig. 2. Bifurcation diagrams: stationary values of the order parameter vs noise self-correlation time at 
e = 0.7 (a = 0.2 (a); a = 0.6, (b)). 

Fig. 3. Phase diagram indicating the ranges of the parameters Tc, e. 

Fig. 4. Basic types of phase portraits in the ordered state: e = 0.7, r = 0.4 (a = 0.2 (a), a = 0.8 (b)). The 
notation for stationary points is the same as in Fig. 2. 

Fig. 5. Dependence of the probability P of realization of various trajectories on the initial value of the order 
parameter xq for e = 0.7 (curves 1, 2 correspond to a = 0.2 r = 0.4, 0.6, curves 3,4 correspond to a = 0.8 
T = 0.4, 0.6. The initial value of the conjugate momentum (j>o = 0.5. 
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